Quantum Charge Liquid - a new metallic state in 2 dimensions ? 
Application to Na .7CoO 2 family 
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'Quantum Charge Liquid'(QCL), a new phase, is proposed to describe anomalous metallic prop- 
erties (T < 100 K) of a density 'y' of doped holes with strong coulomb repulsions, in a band insulator 
NaCo02 with N sites. This phase, a quantum melted charge order, is characterized by a 'charge only' 
fermi-sea (cFS) of yN charges and a 'soft spin-liquid' of yN spin-i moments forming an 'effective 
lattice' of Heisenberg antiferromagnet. An approximate microscopic theory and some consequences 
are outlined. Anomalous 'Fermi surfaces' seen in ARPES of Nao.sCo02 and Nao.7Co02 matches 
well with our cFS, rather than a corresponding single band free electron FS. 
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Metallic solids exhibit a rich variety of symmetry 
broken phases such as, superconductivity, antiferromag- 
netism, ferromagnetism, charge density wave etc. In ad- 
dition they also exhibit a variety of 'normal state phases', 
without any obvious symmetry breaking. A popular and 
ubiquitous normal state in metal physics is Landau's 
fermi liquid state [1] that has been dominating the field 
for the past several decades. With the advent of lower 
dimensional conductors there has been a paradigm shift 
and novel metallic states such as Luttinger liquid in ID 
[2], tomographic Luttinger liquid and a variety of quan- 
tum spin liquids with charges in 2D [3] , quantum critical 
states [4] and others have come to the scene. Theories of 
these new metallic states often involve some exotic quan- 
tum world and modern mathematical ideas; at the same 
time they are believed to have far reaching implications 
for the real world through nano-technology and devices. 

The aim of the present paper is to suggest a normal 
metallic state, in the temperature below about 100 K, 
called quantum charge liquid (QCL) in Nai_ y Co02 , a 
quasi two dimensional metal [5]. Our motivation arises 
primarily, from among a variety of anomalies [6-14], a 
contrasting metal like quasi-particles (with a significantly 
expanded fermi surface !) seen in ARPES and an insu- 
lating type Curie- Weiss spin susceptibility. 

Our working hypothesis is that there is a microscop- 
ically homogeneous single quantum liquid that exhibits 
the above anomalous behavior. We define a QCL state in 
a band insulator containing a density 'y' of doped elec- 
trons (or holes), living in a single band with coulomb 
repulsions. In a QCL, the underlying electrons become 
spin less fcrmions and condense into some kind of charge 
fermi sea (cFS). The yN spin-^ moments 'riding' on the 
charges form a weakly interacting soft (vanishingly small 
spin propagation velocity) non-condensed bose fluid of 
spin-i moments in an effective Heisenberg lattice. Phe- 
nomenologically, QCL in Nai_ y Co02 seem to arise from 
a quantum melting of certain frustrated charge ordered 
states. 

A schematic 'normal state phase diagram' for a 2D nar- 
row band metal is shown in figure 1 , where the suggested 



QCL appears well away from the Mott insulating end. It 
is also interesting to note a kind of duality that exists 
between charges and spins between the ends of the phase 
diagram. In contrast, in a ID narrow band metal we have 
only a Luttinger liquid phase. In the context of cuprates 
there has been attempts, notably by Putikka and collab- 
orators [15] looking for spin less fermions, close to the 
Mott insulating end; they however find hard core boson 
as a more appropriate description, consistent with our 
phase diagram. An earlier work [16], with no long range 
coulomb interaction, found a different normal state phase 
diagram. 
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FIG. 1. Schematic 'Normal State Phase diagram' of a 2D 
narrow band metal. The proposed QCL appears away from 
the Mott insulating end; y is the carrier (holes in the case 
of Nai_ y Co02 ) density. A type of 'duality' between charges 
and spins at both ends is visible. Shadow denotes a crossover. 



Our QCL hypothesis explains the shapes and area of 
the ARPES fermi surfaces [10,12] in a satisfactory fash- 
ion and some of the other features qualitatively, mak- 
ing the idea worth pursuing further. To this end, we 
focus on a one band model introduced by us [17] for 
Na x Co02.yH20 , a new superconductor [14] and formu- 
late the QCL hypothesis first in terms of a variational 
wave function and then a slave fermion representation. 
To gain more understanding and motivate further anal- 
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ysis we also present a Landau type energy functional. 

There are many striking anomalies in the experimental 
results in Nai_ y Co02 for 0.5 > y > 0.25: i) Curie- Weiss 
behavior of uniform spin susceptibility [7-13], and non- 
Korringa ^ behavior [11] in Na-NMR ii) a large ther- 
moelectric power [6], which is easily suppressed [8] by 
small external magnetic fields <~ 10 T hi) temperature 
dependent Hall constant and a high temperature linear 
T behavior [9], iv) p a b ~ Tor x 1_e over a temperature 
range <T < 100 K and low c-axis conductivity [8,9,13], 
v) Large fermi surface radii seen in ARPES [10,12], too 
large compared to the non interacting fermi surface, vi) 
charge orderings as inferred from NMR/NQR [7,11]. 

Our model for Nai_ y Co02 is a single orbital tight 
binding model on a triangular lattice with yN hole carri- 
ers added to a band insulator NaCoC>2. Even though the 
first LDA calculation [18] showed more than one band 
at the fermi level for Nao.sCoC^ , we work with our ear- 
lier single band hypothesis [17], that has support from 
ARPES [10,12] and recent LDA + U calculations [21,20]. 
That is, in our modelling, in NaCo0 2 each Co 3+ site is an 
orbitally non-degenerate spin singlet. The doped 'holes' 
are orbitally non-degenerate Co + sites; they carry a 
spin-i moment and a charge e + , with reference to a Co 3+ 
site. 

Two cases of interest to us are Na .5CoO 2 and 
Na .7CoO 2 , corresponding to hole band at fillings 0.25 
and 0.15. Both cases are away from the Mott insulat- 
ing half filled band. Since the band fillings are low, un- 
screened short range coulomb interactions are important. 
Simple estimates, similar to the one made in reference 
[19] shows that the unscreened interactions are compa- 
rable in strength to the band width, making charge or- 
der very likely. For the commensurate cases y = ^, | 
and j, j we get triangular lattice charge orders: \/3 x \/3 
and 2x2. However, charge order for an arbitrary 'y' 
is frustrated by strong commensurability effects of the 
triangular lattice and an equally important disorder po- 
tential from the dopant Na + layer. We suggest that the 
result is a quantum melting into a novel, nearly uni- 
form and homogeneous metallic state, QCL. The quan- 
tum molten state will be dominated quantum delocalized 
defects (topological and non-topological) of these orders. 
We do not discuss the mechanism of quantum melting, 
but focus on the quantum melted state. 

In our narrow band metal, in addition to the usual su- 
pcrcxchangc, multiparticle exchange processes riding on 
local charge order correlations in the quantum melted 
state are likely to be present. We put them together and 
write an effective single band projected model: 

H t , = -tJ2cj a C ja + H.c. + J]T(S 4 • S, - \n inj ) + 

(ij) {ij) 
ij i 



with a local constraint rij = + ^ 2 of no double 
occupancy of holes; it means projecting out Co 5+ states. 
Here C's are hole operators and S's spin operators. Since 
we use hole operators, in contrast to reference [17], the 
sign of the hopping matrix element is positive. The val- 
ues of various parameters [17] are t m 0.1 eV, J w 7 meV; 
the screened near and next nearest neighbor coulomb re- 
pulsion parameters are Vi w 0.8 eV and V 2 ~ 0.4 eV. 
And e, is a random potential arising from the Na + layer. 
A simple estimate shows that the root mean square fluc- 
tuation of ei can be as large as \J (ef ) w 0.2 eV. In the 
absence of disorder and incommensurability effects the 
above model will exhibit charge order [19] for the range 
of doping < y < 0.5, with an energy scale of the order 
of 100 - 300 K, as the coulomb scale is comparable to the 
band width. Experimentally we see a metallic behavior, 
indicating a quantum molten state; any charge ordering 
signal seen in NMR/NQR [7,11] is likely to be very small 
in amplitude. This needs to be studied further. 

The above is a hard many body problem, containing 
disorder and strong interaction. As a first step in gaining 
some understanding, we propose a variational wave func- 
tion for our phcnomenologically motivated (non-random) 
QCL state in coordinate representation: 

*QCL [ri^] ~ * c h(ri, ..r yN ) Xspin(ricri, ..r yN cr yN ) (2) 

Here r, denotes the lattice site of the hole and cr,, the 
z-component of its spin. The charge function ^ c h is anti- 
symmetric and the spin function x S pin symmetric, mak- 
ing 'J'qcl antisymmetric. We choose \l/ c h to be a simple 
slater determinant of a charge fermi sea(cFS) times a 
Jastrow factor: 

* c h(ri ) ... ) r yN )~n/( r «)E(- 1 ) Npe<E,k, ' PP ' ( 3 ) 

i<j p 

where P stands for permutation of the particle indices, 
Np = 1 (0) for odd (even) permutation. And k; denote 
N y k-points inside the charge fermi sea. We choose a 
form of f(rij) that gives us a quantum melted liquid of 
charges, with some short range charge order correlations. 

The spin part of the wave function is a bose fluid of 
spin-i moments - a soft spin liquid state (RVB) without 
any long range magnetic order 

To discuss the low energy physics of QCL wave func- 
tions we find it convenient to use slave fermion repre- 
sentation, G[ — e-i& ia with anti-commutation relations 
for the fcrmionic charges {e i7 ej} = Sij and commutation 

relation for the bosonic spins [s io ., sj- CT /] = 5ij5 a ^>. Our 
Hamiltonian (equation 1) becomes: 

HtJ - t E e]e iS ls j( , + H.c. - J £ + 

(ij) (ij) 
ij 
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Here — 4g J2 a aC ia C is . An early slave-fermion [23] 
mean-field theory for a 2D t-J model on a square lattice 
(with no short range coulomb interaction) focussed on a 
ferromagnetic solution for nearly all y and J << t. We 
are however interested in understanding anomalous nor- 
mal state and hence do not focus on symmetry broken 
ground states such as charge and or magnetically ordered 
states in the present paper. 

Notice the change in sign of the hopping term in the 
slave-fermion representation, in contrast to the slave- 
boson representation. Within the Hilbert space of the our 
variational wave functions the charge dynamics is easily 
represented by an approximate effective charge (ej,ej) 
Hamiltonian: 

H c h « txo ^ etgj + H.c. 

k v/3 

~ ^[ 2 ^°( cos kx + 2 cos -f C0S ^T k y)~ Ai]4 g fe ( 5 ) 




FIG. 2. Fig 2a) Charge Fermi surface of QCL for different 
fillings (y = 0.5, 0.3, 0.15; green, red and blue). 2b) An ap- 
proximate sketch of ARPES fermi surfaces of reference [10] 
and [12] for y =0.5 and 0.3, which are strikingly close to our 
cFS of figure 2a. 2c) Corresponding free hole fermi surfaces, 
in our single band model (y = 0.5, 0.3, 0.15; green, red and 
blue) . Bottom inset gives density of states and band filling of 
for y = 0.3 case: left - non-interacting hole band with filling 
0.15; right - charge band of QCL with filling 0.3. Notice the 
flip of the charge band and doubling of filled area. 

The Jastrow factor and the ensuing short range charge 
order correlation renormalizes the hopping integral, t — > 
t. The spin-back flow also renormalizes the charge dy- 
namics through xq ~ (Ecr s la s ja)- Both < and Xo may 
be formally computed as many body averages with re- 
spect to the QCL wave function. In the present paper 



we do not attempt to calculate t\o, but consider them 
as parameters to be determined from experiments; we 
estimate this to be i\o ~ 1-10 meV from existing exper- 
iments [8-13]. 

Let us discuss the consequence of the renormalized 
charge fermion Hamiltonian. Because of the sign of the 
hopping term, the band bottom of the charge fermions 
are at the K and K' points in-contrast to the holes whose 
band bottom is at the T point, leading to a flip in the 
density of states (figure 2). The shapes of the charge 
fermi sea (cFS) are shown in figure 2a for three values 
of y = 0.5 (green), 0.3 (red) and 0.15 (blue). It is very 
interesting that our cFS are close to (about 5 to 10 % 
larger) the ARPES results (figure 2b) for y = 0.3 and 0.5 
including the hexagonal distortion that is present for y = 
0.3 and absent for y = 0.5. For comparison the corre- 
sponding free fermi surface of yN holes with spins, in the 
single band model is shown in figure 2c for the three fill- 
ings. If our single band modeling is correct, correlation 
has indeed made a remarkable change. Further the usual 
Luttinger theorem is violated and there seems to be a new 
Luttinger theorem and Luttinger volume containing yN or 
(l-y)N k-points. 

Even though there is a nominal factorization of charge 
and spin subsystems in our wave function (equation 2), 
similar to the Ogata-Shiba wave function [22] in the 
U = +oo ID Hubbard model, actual dynamics is per- 
formed by electrons or holes that carry both charge and 
spin. In view of this, experimentally measured quanti- 
ties or physical (gauge invariant) correlation functions 
are convolutions of certain charge averages and spin av- 
erages, representing the strong back flow constraints ac- 
companying a 'charge' or 'spin' motion. So a nominal 
charge fermi surface that is present in the charge part of 
our variational wave function may be completely washed 
out on convolution. 

In the absence of any hard many body calculation we 
interpret the presence of 'quasi-particle' peak in ARPES 
and its near coincidence with our charge fermi surface, 
to imply that the back flow constraints do not wash the 
charge fermi surface completely. That is, an electron 
Greens function may be approximated by a convolution 
of charge and spin green function. Symbolically, 

G(electron) ss G(charge) x G(spin) (6) 

convolution 

In a spin charge separated system depending on the dom- 
inance of the singularity of G (charge) or G(spin) the 
G (electron) will exhibit the corresponding singularity. 
For example in cuprates, G(spin) has a stronger singular- 
ity than the nearly incoherent and soft G (charge), leading 
to the observed 'electron fermi surface' albeit with power 
law singularities of the spectral functions, reflecting the 
non-fermi liquid behavior. In QCL we expect G(charge) 
to be more singular than the nearly incoherent and soft 



3 



G(spin), leading to the observed 'charge Fermi Surface'. 
It will be interesting to investigate further the non-fermi 
liquid spectral function of electron in QCL using ARPES. 

The spin part of the wave function is a bose fluid of 
spin-i moments - a soft spin liquid state (RVB) without 
any long range magnetic order. Let us discuss the spin 
dynamics within the Hilbcrt space of our variational wave 
functions. Since we have a overdoped case, the superex- 
change is considerably reduced J — > y 2 J (y 2 is the prob- 
ability of finding two spins, Co 4+ on neighboring sites); 
further the nearest neighbor short range repulsion en- 
courages certain multiparticle exchange processes within 
the fluctuating charge ordered regions. The relatively 
fast charge dynamics allows us to average over dominant 
charge configurations and define a fictitious lattice (con- 
taining yN sites) of weakly interacting spin-i moments, 
in the spirit of Born-Oppenhcimer approximation. These 
moments have a Heisenberg antifcrromagnctic coupling 
J c ff «J y 2 + Ji < J. Here Ji represents effective exchange 
couplings arising from other processes. The effective spin 
Hamiltonian of this fictitious lattice containing yN sites 
is : 

-ffspin ~ JcS ^ S; • S TO (7) 

(Im) 

The above spin Hamiltonian is meaningful in getting any 
kind of response of the system at length scales I > > , 

where -7= is the mean inter-hole distance in units of 

Vv 

the lattice parameter. In the same fashion the above 
spin Hamiltonian is also meaningful only for time scales 
t >> the hole derealization time scale. 

As mentioned earlier, for y — \ and i charge ordered 
state are the triangular lattices: \/3 x and 2x2. We 
take these lattices as the lattice of the fictitious Heisen- 
berg spins, in trying to understand fillings close to 0.5 
and 0.7. The above frustrated spin system has interest- 
ing spin liquid physics. 

The high temperature limit of the above spin Hamilto- 
nian may be thought of as a classical gas of spins- i mo- 
ments, suggested in reference [8] that carry the spin en- 
tropy in the thermopower measurements (through back 
flow) and also determines the Curie- Weiss susceptibil- 
ity seen in experiments. In this sense J c ff may be also 
determined from the existing experimental results. Our 
estimate of this from experiments [7-13] is J e ff ~ 1 meV. 

All the above considerations also leads in a natural 
way to a Landau like microscopic energy functional for 
our QCL: 

£[{n k }; {S ; }] w e k h «k + ^ /k,k'n k n k / + 

k 

+ Jeff^(S ; -S ro -i) (8) 

Here e£. h is the energy of the charge carrying quasi- 
particles and f k y are the Landau parameters, which 



may also be singular for the 'chargons'. Various Lan- 
dau parameters can be treated as a phenomenological 
parameters, to be determined experimentally. 

To summarize, the metallic state in Nai- y Co02 with 
reasonably defined quasi particle, as seen in ARPES, and 
coherent transport occurring over a narrow temperature 
interval of about < T < 100 K has been identified with 
a novel Quantum Charge Liquid. Insulating like Curie- 
Weiss spin susceptibility gets naturally explained. We 
also find that the Hall constant in our theory is negative, 
consistent with the experimental observation [9] , over the 
above temperature range. There has been experimental 
signals [24] for magnetic phase transitions at low temper- 
atures in Nai_ y Co02 • QCL can support low tempera- 
ture symmetry broken states, including chiral metallic 
phase [17]. 

We believe that the idea of QCL may have a wider ap- 
plicability: i) certain \ filled 2D organic conductors, ii) 
overdoped cuprates (La2- a; Sr a ;Cu04) with x > 0.4, (if it 
can be synthesized) and iii) metal-insulator transition in 
2D in the presence of interaction and disorder. 

I thank Debanand Sa and R. Shankar for discussions. 
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